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ABSTRACT 

In  this  paper,  an  analysis  is  performed  to  study  the  momentum,  heat  and  mass  transfer  characteristics  of  MHD 
natural  convection  flow  over  a  moving  permeable  surface.  The  surface  is  maintained  at  linear  temperature  and 
concentration  variations.  The  non-linear  coupled  boundary  layer  equations  were  transformed  and  the  resulting  ordinary 
differential  equations  were  solved  by  perturbation  technique.  The  solution  is  found  to  be  dependent  on  several  governing 
parameters,  including  the  magnetic  field  strength  parameter,  Prandtl  number,  Schmidt  number,  buoyancy  ratio  and 
suction/blowing  parameter,  a  parametric  study  of  all  the  governing  parameters  is  carried  out  and  representative  results 
are  illustrated  to  reveal  a  typical  tendency  of  the  solutions.  Numerical  results  for  the  dimensionless  velocity  profiles, 
the  temperature  profiles,  the  concentration  profiles,  the  local  friction  coefficient  and  the  local  Nusselt  number  are 
presented  in  various  combinations  of  parameters. 
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1.  INTRODUCTION 

Hydromagnetic  incompressible  viscous  flow  has  many  important  engineering  applications  such  as  magneto 
fluid  hydrodynamic  power  generator  and  the  cooling  of  reactors.  Tsou  et  al.  [Ill]  Considered  the  effect  of  heat 
transfer  on  a  continuous  moving  surface  with  a  constant  velocity  and  experimentally  confirmed  the  numerical 
results  of  Sakiadis.  Chen  and  Strobel  [IV]  considered  the  effect  of  a  buoyancy  induced  pressure  gradient  in  a 
laminar  boundary  layer  of  a  stretched  sheet  with  constant  velocity  and  temperature.  Chakrabarti  and  Gupta  [V] 
studied  the  MHD  flow  of  Newtonian  fluids  initially  at  rest,  over  a  stretching  sheet  at  a  different  uniform 
temperature.  Vajravelu  and  Hadjinicolaou  [VI]  made  an  analysis  of  flows  and  heat  transfer  characteristics  in  an 
electrically  conducting  fluid  near  an  isothermal  sheet.  Cheng  and  Huang  [VII]  considered  the  problem  of  unsteady 
flows  and  heat  transfer  in  the  laminar  boundary  layer  on  a  linearly  accelerating  surface  with  suction  or  blowing  in 
the  absence  and  presence  of  a  heat  source  or  sink.  Xu  and  Liao  [VIII]  investigated  the  un-steady  MHD  flows  of  a 
non-Newtonian  fluid  over  a  non-impulsively  stretching  flat  sheet  and  presented  an  accurate  series  solution. 
Some  works  are  available  on  the  subject  of  magneto -hydrodynamic  (MHD)  convection  in  porous  medium  was 
studied  by  Abdelkhalek  [IX-XI].  The  problem  of  magneto-hydrodynamic  natural  convection  about  a  vertical 
impermeable  flat  plate  is  presented  by  Wilks  and  Hunt  [XII].  Watanabe  and  Pop  [XIII]  studied  the  case  of  a  wedge. 
Kafoussias  [XIV]  investigated  the  MHD  free  convective  flow  through  a  porous  medium  over  an  isothermal  cone 
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surface.  Bakier  et  al.  [XV]  Presented  non-similar  solutions  for  free  convection  from  a  vertical  permeable  plate  in  porous 
media.  Similarity  solutions  of  natural  convection  boundary  layers  adjacent  to  vertical  and  horizontal  surfaces  in  porous 
media  with  internal  heat  generation  were  studied  by  Pop  [XVI].  Grubka  and  Bobba  [XVII]  investigated  heat  transfer 
characteristics  of  a  continuous,  stretching  surface  with  variable  wall  temperature  Ali  [XVIII]  numerically  presented  the 
heat  transfer  characteristics  of  a  power  law  continuous  stretched  surface  without  and  with  suction  injection.  Takhar  et  al 
[XIX].  Studied  a  MHD  symmetric  flow  past  a  semi  infinite  moving  plate  and  numerically  obtained  the  solutions  Yih  [XX] 
studied  a  free  convection  effect  on  MHD  coupled  heat  and  mass  transfer  of  a  moving  permeable  vertical  surface. 

The  main  objective  of  this  analysis  investigates  the  effect  of  heat  and  mass  transfer  in  a  hydro-magnetic  flow  of  a 
moving  permeable  vertical  surface.  In  the  formulation  of  the  problem,  it  can  be  demonstrated  that  the  system  of 
momentum,  heat  and  mass  conservation  equations  can  be  reduced  to  some  parameter  problem  by  introducing  suitable 
transformation  variable.  The  resulting  transformed  governing  equations  are  then  solved  using  a  perturbation  technique  [I]. 
Typical  results  for  the  velocity,  temperature  and  concentration  profiles  are  presented  for  various  governing  parameters. 
Also  the  local  skin-friction  coefficient  as  well  as  the  heat  and  mass  transfer  results  is  illustrated  for  representative  values  of 
the  major  parameters. 

2.  MATHEMATICAL  FORMULATION 


Let  us  consider  the  two-dimensional  free  convection  effect  on  the  steady  incompressible  laminar  MHD  heat  and 
mass  transfer  characteristics  of  a  linear  moving  permeable  vertical  surface  when  the  velocity  of  the  fluid  far  away  from  the 
plate  is  equal  to  zero.  The  variations  of  surface  temperature  and  concentration  are  linear.  All  the  fluid  properties  are 
assumed  to  be  constant  except  for  the  density  variations  in  the  buoyancy  force  term  of  linear  momentum.  The  magnetic 
Reynolds  number  is  assumed  to  be  small  so  that  the  induced  magnetic  field  is  neglected.  No  electric  field  is  assumed  to 
exist  and  both  viscous  and  magnetic  dissipations  are  neglected.  The  Hall  effect,  the  viscous  dissipation,  and  the  joule 
heating  terms  are  also  neglected.  Under  these  assumptions,  along  with  Boussinesq  approximations,  the  boundary  layer 
equations  for  this  problem  can  be  written  as: 


d u  dv  _ 

dX  +  dX  ~  (2-1) 


du  d2u  0/rr  _  .  0  .  _  _  .  aB2u  uV 

+  =  -TJ  +  g(5c{C-CJ - * - r 

dX  dY  dY  p  a 


dT  dT  d2u 

u - f  v  —  =  an  — - 

dX  dY  0  dY2 


dC'  dC'  „d2C'  \ 

u— —  +  v——  =  D—^—K1(C  -CJ) 
dX  dY  dY2  1 

The  relative  boundary  condition  can  be  written  as 


y  =  0,  v  =  -Vwu  =  BX  ,T  =  Too+aX,C'  =  Coo+bX 
y  — »  °o,  u  =  0,  T  =  ,  C '  = 


(2.2) 


(2.3) 


(2.4) 


(2.5) 
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Where  V  the  uniform  surface  is  mass  flux  and  B,a  and  b  are  prescribed  constants. 


We  introduce  now  the  following  non-dimensional  variables 


dj/_  B  F( 

M  _  >  U  -v  >  ^  \  /  »  *  \tf )  I — 

\v  X^JBv 


0(J1) 


dY  dX 
C'-C' 


v  ,m=-  r“r" 


C„(X)-Cc 


TW(X)-T„ 

-,y/  =  XF(rj)yfBv,u  =  BXF  \v  =  -F(jj)4bv, 


(2.6) 


Ba'  Kx  v  aB 2  pcb  v 

a  = - ,K=-±,SC  =—  ,M  =-±,N  =  tp-,Pr=  — 

v  B  D  Bp  Pa  a0 

With  a  new  set  of  independent  and  dependent  variable  defined  by  equation  (2.6),  Equation  (2.1)  is  identically 
satisfied,  and  the  partial  differential  equations  (2.2)  to  (2.4)  transform  into  ordinary  differential  equations  (2.7)  to  (2.9). 


f 


F '"+  FF 


M  +  —  +  F' 
a 


\ 


F '  =  -GRp  ( 0  +  NC ) 


(2.7) 


0"+P(Fd'-F'0)  =  O 


(2.8) 


C "+  Sc (FC F'C -  KC)  =  0 


(2.9) 


Where  primes  denote  differentiation  with  respect  to  7] .  The  appropriate  flat  plate,  free  convection  boundary 
conditions  also  transform  into  the  form, 

rj  =  0,F  =  FwF'  =  l,e  =  l,C  =  l] 

n^°°,F'=o,e=o,c=o  J  (210) 

The  velocity  components  are 

u=BXF\v  =  -fJb% 

Where  M=  oB20B'1  q'1  is  magnetic  parameter,  GrR  =  g PaB  2  is  the  buoyancy  parameter,  when  GrRc  —  0 

the  governing  equation  is  reduces  to  the  forced  convection  limit.  However,  as  GrRc  — >  °°  free  convection  is  dominated. 

The  buoyancy  ratio  N  measures  the  relation  importance  of  mass  and  thermal  diffusion  in  the  buoyancy  driven  flow.  It  is 
apparent  that  N  is  zero  for  thermal  driven  flow,  infinite  for  mass  driven  flow,  positive  for  thermally  assisting  flow,  negative 

for  thermally  opposing  flow,  Pr  is  the  Prandtl  number  and  Sc  is  the  Schmidt  number,  F  =  V  /  yj Bl)  is  the  suction  / 
blowing  parameter.  In  the  case  of  suction  Vw  >  0  and  Fw  >  0.  In  the  case  of  blowing  Vw  <  0  and  hence  Fw  <  0. 

3.  SOLUTION  OF  THE  PROBLEM 

The  resulting  differential  equation  contains  arbitrary  parameters,  the  Prandtl  number,  the  magnetic  field  strength 
and  buoyant  force,  the  ratio  of  the  Hartmann  number  is  a  measure  of  the  relative  influence  of  the  magnetic  and  buoyant 
forces  on  the  temperature  and  flow  field.  Solution  of  the  resulting  semi-infinite  domain,  non-linear  equations  is 
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accomplished  with  a  three  part  serie’s  method  [II].  The  employed  power  series  equation  (3.1)  contains  a  term  A  that 
satisfies  the  boundary  condition  and  the  differential  equation  at  infinity,  a  second  term  that  satisfied  the  boundary  condition 
at  zero  and  is  the  solution  of  initial  homogeneous  differential  equations  and  additional  terms  that  are  utilized  to  obtain 
increased  numerical  accuracy.  This  accuracy  is  limited  by  the  number  of  terms  that  will  not  initiate  divergence  of  the 
numerical  results. 


F  —  A  +  £FX  3-  £  F2  +  £  F 3  + _ 

6  —  £0X  +  £2Q2  +  £ 3  #3  +  .... 

C  —  £CX  +  £~C2  +  £  C3  + _ 


(3.1) 


Subject  to  the  boundary  conditions 

7  =  0,  Fj(0)  =  Fw,F2(  0)  =  F3(  0)  =  0,  F/(  0)  =  1,  F2(  0)  =  F3(0)  =  o' 

ex  (0)  =  1,  e2  (0)  =  e3  (0)  =  o,  c,  (0)  =  1,  c2  (0)  =  c3  (0)  =  o  > 

7  ->  00,  Fn  (00)  =  0, 6n  (00)  =  0,  Cn  (00)  =  0,  n-  1,2,3 


(3.2) 


Boundary  conditions,  equation  (3.2),  contain  undetermined  parameter  s  which  aides  in  the  collection  of  terms  for 
each  set  of  the  resulting  linear  differential  equations.  In  some  problems,  it  will  have  physical  meaning  which  results  in  a 
power  series  in  that  parameter.  However,  the  present  case  s  equals  unity.  Substitutions  of  the  series  representation  into  the 
differential  equations  and  collection  of  terms  by  like  powers  of  s  result  in  a  family  of  linear  differential  equations,  and  the 
first  three  sets  are 


r 


F{  +  AFX 


1 


M  +  — 

v  aj 


F\  -  ~G,  R~2  (#i  +  NCX ) 


(3.3) 


e\  +  APr0[  =  0  (3.4) 

C,"  +  AScC[  -  KCX  =  0  (3.5) 


f 


F2  +af2 


1 


M  +  — 
a  j 


K  =  K-FA-GrR;2{02  +  NC2) 


(3.6) 


o2+APre2  =  pr(F;el-Fle'l) 


(3.7) 


C’  +  AScC2  - KC2  =  Sc(FxCx  - F{C[) 


F3  +  AF3 


f  O 

M  +  — 

V  CCj 


K  =  2FF  -  f2f;  g,r (e, + 1 vc, ) 


e\  +  APre,  =  pr(F:e,  +  F,e,  -F&-  F,e[) 


C3  +  ASC3  —  KC,  —  S c ( F| C2  +  F2Cx  —  FxC2  ~  F2CX) 


(3.8) 

(3.9) 

(3.10) 

(3.11) 
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On  solving  equations  (3.3)  to  (3.11)  by  using  boundary  conditions  (3.2),  we  get  the  following  solutions 


ex  =  e-AP'v 

(3.12) 

C,  =  e~aiT1 

(3.13) 

Fj  =  a5  +  a4e~a 3n  +  a6e^APrV  +  a7e~ain 

(3.14) 

02  =  ane-AP'n  +  a^APM  +  a9eH,h+AP 

(3.15) 

C2  =  a13e“w  +  aneHa2+a3)T1  +  aX3e<a2+AP')n  + 

(3.16) 

F,  =  a36e~a3'1  +  a37  +  a24e~(a3+APr)r’  +  a25e~<'a2+APr)>1  +  a26eHa2+“3)n  +  a27Tje~c‘3'1 
+ai4e~c‘1'1  +  a29rie~airi  +  a35e~APr’1  +  a32rje~APr,} 

(3.17) 

n  -AP./7  .  -(a?+APr)n  ,  __  ~(a-i+APr)n  .  -(a7+APr)rj  .  „  -(a-y+APr)n 

03  =  a64e  rl +a5le  3  rl +a57rje  3  r,+a53e  2  T‘  +a54rje  2  r/ 

.  -(dn+an+AP,.)?!  ,  -(2 a?+APr)Tj  .  -(a-,+2APr)rj  .  -(2a7+APr)ri  ,  -(a7+2AP.)ri 

+a55e  2  3  rl+a56e  3  r  l+a51e  3  r  l+a58e  1  r  l+a59e  2  r  l 

i  -2APn  .  „  -2APrn  .  „  ~APrn  .  „2  -APn 

+a60e  1  +abir/e  ''+a62Tje  ,l+a63rje 

(3.18) 

r  -n  p~a2jl +a  p~(ai+a3)'1 +a  n +  a  P<a^a^APCv  +n 

+aue~(a2+AP')v  +a,2T]e-(a2+AP')’1  +  a,3e~{a2+2AP^  +  aMe~(2a2+a^  +  aS5eH2a2+AP^ 

+a36e~2a2n  +  a37r]e~a2n  +  a8S//2  e~ain 

(3.19) 

r-<  ~<hTI  i  i  -(a^+APr)7j  .  „  -(a-,+APr)rj  .  -(a>+APr)n  .  „  -(a7+APr)n 

F3=al42e  31  +aU3  +  an6e  3  rl +aul7je  3  r  ,+ame  2  r  / +an9r/e  2  r  ' 

.  -(a?+ao+AP)n  .  -(a7+a?)7j  .  „  -(a7+a-,)ri  .  -26077  ,  -260/7  .  -26077 

+ai2Qe  -  +  ame  +al22rje  2  3  l+ame  31  +an4rje  31  +a]25e  21 

.  -(2 6o+APr)77  .  -(6o+2APr)77  .  -(2 a->+APr)n  .  -(a7+2AP.)n  .  -(cp,+2a,)n 

+al26e  3  r,/+al27e  3  r,+ame  2  rl +al29e  2  r'+ame  -  31 

+ame~2APr 11  +  amr/e^2APrr>  +  al33e~t2a2+a3),>  +  al34e~AP'n  +  al33r/e~APr'’  +  al36rj2e~AP,r’ 

+ax37e~a2ri  +  an&T]e~a2ri  +  al39Jj2e~aiV  +  auj]e~a3n  +  al4lij2e~“3’1 

(3.20) 

9[  =  - APreAP » 

(3.21) 

C,  =  -a^l,2n 

(3.22) 

Fx  =  -a3a4e^a3V  -  APra6e^AP,’]  -  a2a7e^“2V 

(3.23) 

02  =-APraue~AP>n  -( a3  +  APr)ase~(a3+APr)'’  ~{a2  +  APr)a9e~('a2+APr)n 
+awe-AP«  -  APrawTje~APr’1 

(3.24) 

C2  = -a2al5e~a2T'  ~{a2  +  a3)al2e~ia2+a3),>  -( a 2  +  APr)al3e~('a2+APr)>1 
+cix4e~ai11  -  a2ax4rje~a2ri 

(3.25) 

www.tjprc.ors 


editor@tjprc.  org 


16 


Uma  Shanker,  Hemant  Poonia  &  S.  S.  Dhayal 


F2  =  -a3a36e  a 3n  -(a3  +  APr)ci24e  ia> +APr)n  ~(a2  +  APr)a25e  <fl2+AP'-)7 


-(a2  +  a3)a26e 
-<r/2a29/7e“fl2'7 


<¥7 

-(02  +a3  >7 


A^,a35e 


*25  C 

+  a27e“'¥7  -  a2a217je~a3n  -  a2a34r/e~a2J1  +  a29e~aiJ1 

-APrrj  .  ~APrlj  a  r>  -AP..1J 

T '  +  a32e  r  -  APra32JJe 


(3.26) 


6»3  =  -APraMe 


-APrlJ 


(a3  +  A/^,)a51e 

^53' 


-(.aj+AP,)!)  _|_  ^  e~(a3+APr)T] 


(a3  +  APr ) 

a52 rje~('cl3+APr)tl  - (a2  +  APr)a53e~(a2+APr)v  +  a54e^a2+AP')ri  - ( a :2  +  APr)a54 

~(a2+APr )//  _  (cii  +  ^  +  APr)a^e-(^a^APr)n  _  +  Ap)a^(2o3+APr)V 

-(a,+2APr)V  _n„  xJPU  iO-(202+iU>r)7_(a^+2Ap)a 


7/e 


-(a3  +  2 Ai^  )a57e  ta3+z/1,;  1,7  -  (2a2  +  AP  )a58e 


~(a2+2APr)j] 


59 


o  A  r>  -2APrn  .  -2AP.n 

■2APra60e  *  +  a6le  ■ 


-APra62T]e  r  7  +  2 ab3rje  r  1  -  APra63rj  e 

C3  =  -a7awe~a2<1  -(a,  +  a3)a77e~<fl2+a3)?7 


'  2APra6lJ]e~2APrV  +  a62e~APrl? 

2-APrn 


-(a2+a3)r/ 


+  (/7g(’ 


-(a2+a3)t] 


~(a2  +  a3)a 


lie 

-(a,  +APr)aue 


-( a2  +c/3  +  APr)a19e  ^-+a^+AP')’1  ~(a2  +  2 a3)a80e 


2  1  “3 

( a2+APr)t i 


+  6(g2£ 


~(a2+APr)tj 


78 

-(a2+2a3)/7 

ta2+APr)j/ 


-(a2  +2 APr)a&3e  (-ai+2ApP,]  -(2 a2  +  a3)ag4e 


a^e~(2a'l+AP,)'1  -  2 a2a86e“2fl2'7  +  a^"*2'7  -  a2a^e~a 2”  +  2agg/7£> 


-(a2  +  APr)a&2T]e 

-{la2+a3)n  -(2ci2  +  APr) 
«2  n 


-a2assrj  e 


2  yJ~aiV 


(3.27) 


(3.28) 


F3  =  -a3a142e 

aulT]e-(a3+AP')T> 

r/e 


-a2n 


(a3  +  APr)ame 


-(a, +APr)t] 


+  ((117^ 


~(a3+APr)t] 


■( a3  +  APr ) 


’  (^2  1  ‘•‘119 

-(a2+a3+APr)^ 


, -(02^)7  ,^-(02^)7  _(  APr)fl 


119 


(a2+APr)n  +a ^  +  APr)a120e  -(a,  +o3)<3121e 


~(a2+a3)r/ 


+al22e 


-(a2+fl3)l? 


-(a2  +a3)al22r/e  <a2+“3)'7  -2a3a123e  2"3'7  +a124e  ^377  -2c/3a124 


2a37]  —  O/y  ^2  /? 

^**2^125^ 


2 

-2a2T] 


(2a3  +  APr)ame 


~(2a3+APr)Tj 


(a3  +  2  APr)avle 


-(a3+2APr)i] 


-(2  a2  +AP)a,7„e 


-(2  a2+APr)rj 


-2 APral3le  ^APrn  +  al32e  2AP,V 


~(a2+2APr)al29e 

-2APrJj 


~(a2+2APr)rj 


~(a2  +  2  a3)a130e 


-(a2+2a3)r/ 
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The  series  for  0,  its  first  derivative  0  TO)  (the  wall  temperature  gradient),  F  '  (the  velocity  profile)  and  F  "(0) 
(the  wall  velocity  gradient)  were  evaluated  and  are  presented  in  Figures  1-7.  The  general  results  of  the  investigations  are 
the  imposed  magnetic  field  diminished  the  velocity  field;  wall  shear,  flow  rate,  and  wall  heat  transfer,  also  the  onset  of  free 
convection  are  retarded  while  the  fluid  temperatures  and  the  time  required  for  the  flow  to  each  steady  state  are  increased. 
In  addition,  sizable  influences  on  the  flow  and  thermal  fields  can  be  produced  with  moderate  magnetic  field  strengths  only 
for  liquid  metal  flows  while  the  effects  of  induced  magnetic  fields  and  Joule  heating  are  very  small.  The  magnetic  field 
strength  is  to  reduce  the  values  of  wall  shear  stress  regardless  of  blowing  and  suction  strength.The  magnetic  force  is  known 
to  have  a  retarding  effect;  it  decreases  the  wall  shear  stress.  Blowing  has  a  similar  effect  which  used  to  aid  the  magnetic 
field  in  reducing  the  wall  shear  stress.  However,  section  has  an  opposite  effect  and  increases  the  shear  stress  at  the  wall 
opposing  the  magnetic  effect.  In  case  of  free  convection,  suction  decreases  wall  shear  stress,  this  may  be  because 
increasing  suction  will  decrease  the  flow  velocity  and  in  turn  diminishing  the  retarding  effect  of  the  magnetic  field  force. 
In  order  to  get  the  physical  insight  into  the  problem,  numerical  calculations  are  carried  out  for  different  values  of  Prandtl 
number,  Hartmann  number,  Schmidt  number  and  Sherwood  number. 
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Figure  1:  Effects  of  Sc  and  a  on  Tangential  Velocity  Profiles,  taking 
Pr=0.72,  Re=50,  M=l,  K=5,  N=20,  Gr=5,  A=l,  Fw  =0,  £  =1 
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Figure  2:  Effects  of  M  and  N  on  Tangential  Velocity  Profiles,  taking 
Pr=0.72,  Sc=0.2,  Re=50,  K=5,  a  =1,  Gr=5,  A=l,  Fw  =0,  £=1 

Figure  1  illustrates  the  influence  of  Schmidt  number  Sc  and  Porosity  parameter  a  on  the  velocity  profile  in  the 
boundary  layer  against  T] .  It  shows  that  the  velocity  profile  decreases  exponentially  as  increasing  t] .  It  is  also  observed 
that  the  velocity  profile  decreases  as  increasing  Schmidt  number.  This  causes  the  concentration  buoyancy  effects  to 
decrease  yielding  a  reduction  in  the  fluid  velocity.  The  velocity  profile  also  decreases  as  increasing  Porosity  parameter. 
The  presence  of  a  porous  medium  increases  the  resistance  to  flow  resulting  in  a  decrease  in  the  flow  velocity. 
This  behavior  is  depicted  by  the  decrease  in  the  velocity  as  a  decreases  as  shown  in  figure  1. 

Figure  2,  illustrates  the  influence  of  the  Hartmann  number  M  on  the  velocity  profiles  in  the  boundary  layer. 
Application  of  magnetic  field  of  an  electrically  conducting  fluid  gives  rise  to  a  resistive  type  force  called  the  Lorentz 
force.  This  force  has  the  tendency  to  slow  down  the  motion  of  the  fluid  in  the  boundary  layer.  Also,  the  effects  on  the  flow 
and  thermal  fields  become  more  so  as  the  strength  of  the  magnetic  field  increases.  This  is  obvious  from  the  decreases  in  the 
velocity  profiles  presented  in  Figure  2.  It  is  also  seen  from  Figure  2  that,  as  expected,  the  velocity  profiles  F '  increase 
monotonically  with  increasing  the  buoyancy  ratio.  The  decrease  of  F '  happens  because  of  an  accelerating  force, 
which  acts  in  a  direction  parallel  to  the  T] -axis.  As  shown,  the  velocity  increases  near  the  plate  with  increasing  bouncy 
ratio.  This  is  because  the  effect  of  the  bouncy  ratio  N  is  to  increase  the  surface  heat  and  mass  transfer  rates. 


Figure  3:  Effects  of  N  on  Non-Dimensional  Surface  Velocity  Gradient,  taking 
Pr=0.72,  Sc=0.2,  Re=50,  M=l,  K=5,  N=20,  a  =1,  Gr=5,  A=l,  F  =0,  £  =1 
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Figure  4:  Effects  of  Pr  and  Fw  on  Temperature  Profiles,  taking 
Sc=0.2,  Re=50,  M=l,  K=5,  N=20,  a  =1,  Gr=5,  A=l,  Fw  =0,  £  =1 

The  effects  of  buoyancy  ratio  N  on  transport  properties  are  illustrated  in  Figure  3.  The  skin  friction  coefficients 
reduce  with  increasing  N.  The  surface  shear  stress  increases  with  the  magnetic  parameter  M  due  to  enhanced  Lorentz  force 
which  imparts  additional  momentum  into  the  boundary  layer.  This  reduces  the  boundary  layer  thickness  which  in  turn, 
increases  the  surface  shear  stress. 

Figure  4  shows  the  temperature  Q  (jj  )  profiles  across  the  boundary  layers  and  different  values  of  the 

suction/injection  parameter  F  and  the  Prandtl  number  Pr.  As  mentioned  suction  corresponding  to  F  >  0 ,  injection  to 

F  <  0  and  F  =  0  to  impermeable  plate.  It  is  known  that  the  imposition  of  wall  fluid  suction  reduces  both  the 

hydrodynamic  and  thermal  boundary  layers  which  indicate  a  reduction  in  the  temperature  profiles.  However,  the  exact 
opposite  behavior  is  produced  by  the  imposition  of  wall  fluid  blowing  or  injection.  These  behaviors  are  clear  from  Figure 
4.  It  is  also  observed  that  the  magnitude  of  temperature  for  air  (Pr=0.72)  is  greater  than  that  of  water  (Pr=7.0).  This  is  due 
to  fact  that  the  thermal  conductivity  of  fluid  decreases  with  increasing  Pr,  result  a  decrease  in  thermal  boundary  layer 
thickness. 

Figure  5  presents  the  non-dimensional  wall  temperature  gradient  #'(0)  local  Nusselt  number  —0 '( M .  0)  at 
selected  values  of  N.  Increasing  the  buoyancy  ratio  will  cause  the  enhancement  in  the  heat  transfer  at  small  values  of  M. 
Increasing  N  will  cause  the  non-dimensional  wall  temperature  gradient  0  '(0)  .  The  local  Nusselt  number  is  found  to  be 
decreased  more  rapidly  at  a  large  value  of  N.  Figure  5  shows  clearly  that  the  presence  of  the  magnetic  field  gives  a 
reduction  of  the  local  Nusselt  number,  in  terms  of  6  \M  ,0)  In  a  manner  that  the  local  Nusselt  number  is  to  be  decreased 

as  M  increases.  Figure  5  illustrates  the  effect  of  mass  transfer  F  on  the  dimensionless  temperature  gradient  profiles. 
It  should  be  noted  here  that  positive  values  of  Fw  indicate  fluid  suction  at  the  surface  while  negative  values  of  Fw 

correspond  to  fluid  blowing  or  injection  at  the  wall.  The  effect  of  suction  is  to  make  the  velocity  and  temperature 
distribution  more  uniform  within  the  boundary  layer.  Imposition  of  fluid  suction  at  the  surface  has  a  tendency  to  reduce 
both  the  hydrodynamic  and  thermal  thickness  of  the  boundary  layer  where  viscous  effects  dominate.  This  result  in 
increasing  the  non-dimensional  surface  temperature  gradient  is  shown  in  Figure  5.  On  the  other  hand,  the  opposite  result  is 
obtained  from  fluid  blowing  Fw  <  0  at  the  wall,  a  decrease  in  the  non-dimensional  surface  temperature  gradient. 


www.tjprc.ors 


editor@tjprc.  org 


20 


Uma  Shanker,  Hemant  Poonia  &  S.  S.  Dhayal 


N,  Fw 
^^-0,  0 

—a— 300,0 
—a— 500,0 
— 0,  -1 
— 0,  -0.5 
-o—O,  0.5 

— t — 0,  1 


Figure  5:  Effects  of  N  and  Fw  on  Non-Dimensional  Wall  Temp.  Gradient, 
Taking  Sc=0.2,  Re=50,  M=l,  K=5,  a  =1,  Gr=5,  A=l,  £  =1,  Pr=0.72 


Figure  6:  Effects  of  K  and  Sc  on  Concentration  Profiles,  Taking 
Pr=0.72,  Re=50,  M=l,  N=20,  a  =1,  Gr=5,  A=l,  Fw  =0,  £  =1 


Figure  7:  Effects  of  Sc  on  Sherwood  Number,  Taking  Pr=0.72, 

Re=50,  K=5,  N=20,  CC  =1,  Gr=5,  A=l,  Fw  =0,  £  =1 

Figure  6  illustrates  the  influence  of  Schmidt  number  Sc  on  the  concentration.  As  Schmidt  number  Sc  increases, 
the  mass  transfer  rate  increases.  Hence,  the  concentration  decreases  with  increasing  Sc.  It  is  evident  from  this  figure  that 
the  concentration  C  takes  its  limiting  value  C*,  for  higher  values  of  the  dimensionless  distance  T|.  From  this  figure, 
we  observe  that  when  the  concentration  difference  SC  is  maintained  constant,  the  dimensionless  concentration  profile 
decreases,  in  the  sense  that  the  values  of  the  Schmidt  number  increases.  It  is  also  concluded  that  concentration  profile 
decreases  as  increasing  chemical  reaction  parameter  K. 
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The  effect  of  Schmidt  number  on  the  Sherwood  number  is  depicted  in  Figure  7.  It  is  clearly  seen  that  the  influence 

of  Sc  of  the  local  Sherwood  number  is  considerable.  The  results  of  figure  indicate  that,  a  larger  SH  is  experienced  for  a 

system  with  a  larger  Sc.  That  is,  the  mass  transfer  rate  increases  with  the  increase  in  Schmidt  number.  This  is  due  to  the 

fact  that  a  larger  Schmidt  number  corresponds  to  a  small  binary  diffusion  coefficient  for  a  given  mixture  and  to  a  thinner 

concentration  boundary  layer  thickness  relative  to  the  flow  boundary  layer  thickness.  This  results  in  a  larger  mass  transfer 

rate  at  the  duct  wall  or  a  larger  Sherwood  number. 

5.  NOMENCLATURE 

T0  -  The  temperature  of  the  plate,  K  ■  T  -  the  temperature  in  the  boundary,  K 
Tw  -  The  temperature  at  wall,  K  ■  Tx  -  Free  stream  temperature,  K 
U,  V  -  Velocity  components  in  X  and  Y  direction  respectively,  m.S  1 
C'  -  Concentration  in  the  fluid,  mol. ill  3.  C'  -  Concentration  at  the  wall,  mol. ill 
-  Free  streaming  concentration,  mol. ill  3 .  C  -  Dimensionless  concentration 
t'  -  The  time,  .S'  .  B,  a,  b  -  constants;  g  -  the  acceleration  due  to  the  gravity,  m.S  1 
Cp  -  The  heat  capacity  of  the  fluid,  J  .Kg  1  .K  .  B(]  -  the  magnetic  induction 


a0  = - (Thermal  Diffusivity);  P  -  the  Prandtl  number; 

pcp 

K,  -  The  chemical  reaction  parameter;  K  -  Non-dimensional  chemical  reaction  parameter 
M  -  The  magnetic  parameter;  Re  -  the  Reynold  number;  N  -  the  buoyancy  ratio 
G  -  The  thermal  Grashof  number;  Sh  -  the  Sherwood  number;  Sc  -  the  Schimdt  number 

2  —1 

D  -  Mass  diffusion  coefficient,  m  .S  ■  F  -  the  velocity 

Greek  Letters 

I  ’  -  The  porosity  parameter;  □  -  non-dimensional  porosity  parameter 
/ 3  -  The  volumetric  coefficient  of  thermal  expansion,  K  1 
f 3  -  Volumetric  coefficient  of  expansion  with  concentration,  K  1 

p  -  The  density  of  the  fluid.  Kg. ill  9  -  Dimensionless  temperature 

2  —1 

p  -  The  coefficient  of  viscosity,  Pa.s;  V-  the  kinematics  viscosity,  m  .S 
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(7  -  The  electrical  conductivity;  K  -  the  coefficient  of  thermal  conductivity.  Win  1 K  1 

6.  CONCLUSIONS 

In  this  paper,  the  effect  of  heat  and  mass  transfer  in  a  hydro-magnetic  flow  of  a  moving  permeable  vertical  surface 
was  analyzed.  In  the  formulation  of  the  problem,  the  system  of  momentum,  heat  and  mass  conservation  equations  were 
reduced  to  some  parameter  problem  by  introducing  suitable  transformation  variable.  The  resulting  transformed  governing 
equations  were  then  solved  using  a  perturbation  technique.  Results  on  the  different  profiles  were  presented  for  various 
governing  parameters.  We  observed  that  the  presence  of  a  porous  medium  increases  the  resistance  to  flow  resulting  in  a 
decrease  in  the  flow  velocity  and  the  velocity  profiles  F '  increase  monotonically  with  increasing  the  buoyancy  ratio  due 
to  increase  the  surface  heat  and  mass  transfer  rates  and  decrease  in  the  increase  of  the  Hartmann  number  M. 
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